Magnetization Plateaus of SrCu2 (303)2 from a Chern-Simons Theory 
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The antiferromagnetic Heisenberg model on the frustrated Shastry-Sutherland lattice is studied 
by a mapping onto spinless fermions carrying one quantum of statistical flux. Using a mean-fleld 
approximation these fermions populate the bands of a generalized Hofstadter problem. Their filling 
leads to the magnetization curve. For SrCu2(B03)2 we reproduce plateaus at 1/3 and 1/4 of the 
saturation moment and predict a new one at 1/2. Gaussian fluctuations of the gauge fleld are shown 
to be massive at these plateau values. 



Two-dimensional (2D) quantum spin systems that do 
not order magnetically at zero temperature are currently 
a subject of great theoretical and experimental interest. 
The recently discovered compound SrCu2 (603)2 is an 
antiferromagnet (AF) with localized spins S* = ^, a gap 
above the ground state ||^ and the unique property that 
its magnetization curve has plateaus at 1/3, 1/4, and 
1/8 of the full saturation moment ||^. The spin system 
may be described by a 2D Heisenberg model on a square 
lattice with exchange constant J' and additional diagonal 
bonds J on half of the square plaquettes (see inset of 
Fig. 0). 

This lattice was studied many years ago by Shastry and 
Sutherland [Q who noted that there is an exact eigen- 
state which is obtained by putting singlets on all diag- 
onal J bonds. This eigenstate is the ground state for 
a wide interval of J'/J. For J'/J smaller than ^ 0.7 
1^ , the system has dimer long-range order and for larger 
J' it has conventional Neel long-range AF order. There 
may be additional phases in between such as a plaquette 
singlet phase |^-^ but they are apparently not realized 
in SrCu2(B03)2 where J'/J is estimated to be smaller 
than 0.65 |^,|o|. This dimer ground state explains the 
spin gap as seen in experiments. However the existence 
of plateaus has no immediate explanation since the sim- 
plicity of the ground state does not extend to the excited 
states of this peculiar lattice. In this Letter we use a 
mean-field approximation to a Chern-Simons (CS) field- 
theoretic approach to quantum magnets suggested some 
time ago ]lT| , p2[ and we obtain an excellent quantitative 
fit of the magnetization curve (see Fig. (|l|)) for realistic 
values of the exchange constants. As we discuss, this may 
be evidence for unconventional character of the plateau 
ground states. 

Starting from the pure dimer state, the first excited 
state can be constructed by first breaking a singlet bond 
into a triplet state. Due to the exchange J', such a state 
will move and we expect a dispersive band of triplets as 



low- lying states at least for strong J'. 
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FIG. 1. Comparison between the magnetization curve of 
SrGu2 (603)2 measured by Onizuka et al. (dashed line) and 
the mean-field result (solid line). Inset: Shastry-Sutherland 
lattice. The exchange interaction is J on black links and J' 
on the dotted ones. 

However due to the peculiar triangular coupling be- 
tween the diagonal bonds, the hopping of the triplet is 
forbidden at low orders in perturbation theory As a 
consequence, the triplet band is very flat, a striking fact 
observed by neutron scattering experiments jT^ . Since 
the triplets are very massive particles, it is natural to ex- 
pect that they can crystallize at finite density, and it has 
been proposed that the plateaus are Wigner crystals of 
triplets [p^|-p^. There exists a spin model which is de- 
rived from the Shastry-Sutherland Hamiltonian p^ ] for 
which the plateaus are demonstrated to originate from 
such ordered states. However in this model there are 
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plateaus at 1/4, 1/2 and 3/4 and the overall shape of 
the magnetization curve is not in agreement with exper- 
iments. A closely related physical picture is obtained by 
describing the magnetized triplets by hard-core bosons 
| p5[ . Then the repulsion may favor charge-density wave 
states that are amongst known insulating phases of the 
lattice Bose gas. 

We take a different approach here, mapping the spin 
problem onto a hard-core boson problem and then solv- 
ing the hard core constraint exactly by a further map- 
ping onto spinless fermions coupled to a CS gauge field 
p^ , p^j20) . Within a mean-field approximation the spin 
excitation gap that produces the observed magnetization 
plateaus arises from some of the Landau level gaps in the 
integer quantum Hall effect for the fermions on a lattice. 

The Hamiltonian for the AF Heisenberg model on the 
Shastry-Sutherland lattice is: 



H Jij Si ■ Sj - SI 



(1) 



where St are spin-^ operators, the exchange couplings 
Jij are equal to J' when i,j are nearest-neighbors on 
the square lattice and equal to J when i, j are related by 
a diagonal bond (J. J' > 0) and the external magnetic 
field B is applied along the z-axis. We then map the spin 
operators to hard-core boson operators : 



H — Hxy + Hz , 



{ij-> 



2tt 
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It is important to note that the real magnetic field 
B applied to the spins acts as a chemical potential for 
the fermions as seen from Eq. (Q) but does not con- 
tribute to the statistical flux. We treat the Ising term 
in Eq.(Q) by a simple mean-field decoupling so it be- 
comes a simple function of the magnetization. The ki- 
netic energy term H^y leads to a Hofstadter ||2l[] problem 
for fermions hopping on the Shastry-Sutherland lattice. 



We use a flux attachment choice leading to flux 0/2 on 
each triangular plaquette. The one-body problem from 
Hxy can be straightforwardly analyzed for rational values 
(f> = 2TT{p/q). 

Fixing the magnetization M gives us the flux and the 
number of fermions through Eq.(H). For this value of 
(j) we compute the band spectrum of H^y and fill the 
bands with the available fermions. The energy of the 
filled bands leads to a first contribution Exy{M). We 
then add the contribution from the Ising interaction Hz 
to obtain the total energy E{M). The magnetization is 
obtained as a function of B by minimizing E{M) — BM. 
The Hofstadter diagram for J = J' is given in Fig. || 
where the lowest curve marks the Fermi level (highest 
occupied state) and the upper one marks the lowest un- 
occupied level. Jumps of the Fermi energy as a function 
of M leads to discontinuity of the slope of the function 
E{M). These jumps corresponds to plateaus in the mag- 
netization curve. The structure of the Hofstadter butter- 
fiy thus reflects itself in the appearance of plateaus. The 
effects of geometry are encoded in the Hofstadter spec- 
trum, whose complexity arises from the diffraction of the 
cyclotron orbits by the lattice. 



(2) 
(3) 



= E '^^M - V2)K- - 1/2) - (n,; - 1/2) , (4) 

where = Sf + 1/2 is the occupation number of site i. 
The bosons are then treated as fermions with an attached 
flux tube carrying one flux quantum of flctitious mag- 
netic field. This can be formulated as an exact mapping 
between the spin problem and spinless fermions inter- 
acting with a statistical CS gauge fleld. In a mean-fleld 
treatment the flux tubes are smeared out into a uniform 
background magnetic fleld. The flux per square plaque- 
tte 4) is then tied to the density of fermions and thus to 
the magnetization of the spin system : 
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0.2 0.3 0.4 0.5 

*/(27t) 



FIG. 2. Hofstadter spectrum for the Shastry-Sutherland 
lattice at J = J' = 1. Vertical lines mark the energy bands 
as a function of the statistical flux per square plaquette. 
Hall conductances cTxy (TKNN integers) are indicated for the 
regions of the spectrum which are explored in the magneti- 
zation process from M/Msat = -1 ((^ = 0) to M/Msat = 
(0 = ^). 

When J is set to zero, the model reduces to the square 
lattice which has Neel long-range order. The wavefunc- 
tions obtained in the spatially uniform CS mean-fleld 
approximation certainly do not encompass this physics. 
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However the magnetization curve obtained from this ap- 
proach is quahtatively similar to that of the ordered sys- 
tem : it is featureless all the way to full saturation (see 
the J = curve in Fig. If we consider the case of 
the triangular lattice, the uniform mean-field leads to 
curve (a) in Fig. ^. In addition to the zero-field gap 
found by Yang, Warman and Girvin [ p^ , there are many 
plateaus in (a). This is unrealistic : the triangular lat- 
tice spin system is known to be long-range ordered and 
hence gapless and its magnetization process shows only 
a single plateau at M/Msat = 1/3. There is a way 
to improve these results by allowing the mean-field to 
have a three-sublattice structure: one introduces three 
fermion densities ua , , nc and numerically searches for 
a self-consistent solution where the flux 4>a matches the 
density Ua on each sublattice. For M = we find that 
the self-consistent solution remains uniform. However, 
for non-zero magnetization the translation symmetry is 
broken. This non-uniform mean-field solution leads to 
a magnetization curve shown in Fig. ^ which is much 
closer to the truth albeit the zero-field ground state re- 
mains unrealistic (Neel long-ranged order is absent) . The 
1 /3 plateau has a semiclassical origin : It is easily under- 
stood in the Ising limit >> J limit and survives up 
to the isotropic point = J |^^. Our calculation re- 
produces these features and we find the three sublattice 
magnetizations ua = ns = 0.922 and nc — 0.155. 




FIG. 3. Magnetization curves for the J ~ J' Shas- 
try-Sutherland model (uniform mean-field). From left to right 
J = (dashed) , 0.75 (fuU), 1.5 (dashed), 2.5 (full), 3.5 
(dashed) and 5 (full). All these results are obtained for J' — 1. 

We have explored the magnetization process of the 
Shastry-Sutherland lattice as a function of the ratio J/ J'. 
The curves M[B) are drawn in Fig. ^. For small J we 
are close to the square lattice result, i.e. a smooth curve. 
Increasing J we observe plateaus developing, the more 
complex structure being in the regime J/ J' w 2 — 3. 




FIG. 4. Magnetization curve of the triangular an- 
tiferromagnet. Full line: mean-field approximation with 
three-sublattices. Dashed line: uniform mean-field. 

For very large J the curves are again simple. There 
remain only two plateaus at and 1/2. The plateau at 
zero field is due to the fact that for large J the tight- 
binding bands separates into two groups with very small 
dispersion and there is essentially a huge gap in the Hof- 
stadter spectrum that does not depend much upon flux. 
For J' = and an appropriate choice of gauge our wave- 
function exactly reproduces the dimer limit. In the inter- 
mediate regime, the curve depends on the details of the 
spectrum. The plateaus have a finite domain of stability 
which is given in Fig. ^ 




FIG. 5. Width of the magnetization plateaus as a function 
of J / J' for M/Mssx — 0, i, i and \. Additional plateaus at 
fractions ^ for n > 5 also exist in the vicinity of J / J' ~ 1.5. 

To reproduce the qualitative shape of the experimen- 
tal magnetization curve for SrCu2 (603)2 we find that 
J' = 29. 5K and J = 74K, the resulting fit is shown in 
Fig.0. While the zero-field gap is not well reproduced, we 
find 13K instead of 34K, the field strengths at which the 
plateaus occur and the roundings close to the plateaus are 
in very good agreement with the experiment. The values 
of J, J' are reasonably close to the recent estimates pO| ] 
from neutron scattering data J' = 43K and J = 71K. 
These values are not precisely known since there is some 
amount of three-dimensional dispersion of the low-lying 
triplet that should be taken into account. We predict 
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prominent plateaus at 1/3 and 1/2 and nothing else un- 
til full saturation. Note also that at 1/4 there is in fact 
an avoided plateau : the value of J' I J is just outside 
the range of stability of the 1 /4 plateau on Fig. I The 
result of CS mean-field theory is quite different from the 
semiclassical analysis of the effective Bose gas |l^ and 
from the exactly soluble spin model [ pTf . 

As can be seen in Fig. || the spin gap (plateau at 
M/Msat — 0) opens at J/ J' = in this approach instead 
of opening at the correct critical coupling J /J' ~ 1.5. 
This feature comes from the fact that the Neel state of the 
square-lattice antiferromagnet is not correctly described 
in the uniform mean-field approximation, as mentioned 
above in the case of the triangular-lattice antiferromag- 
net. This is corrected by computing non- uniform solu- 
tions of the mean-field CS approach with two-sublattices 



From the good fit in Fig. |] we deduce that the CS 
mean-field theory works better for magnetized than un- 
magnetized states. This is because the initial mapping 
to hard core bosons selects out a preferred spin quantiza- 
tion axis which, together with the mean-field treatment 
of the Ising term, obscures the SU(2) symmetry in the 
underlying spin problem. The Zeeman energy reduces 
the SU(2) to the U(l) symmetry associated with conser- 
vation of S^, which in the language of the fermions is 
simply particle number conservation. We note that the 
physical magnetic field breaks T symmetry in the orbital 
part of a Hubbard model description of the system but 
has no effect in a pure spin model description other than 
to introduce a Zeeman term. However our mean field so- 
lution does break T symmetry. The increased size of the 
'magnetic' unit cell (due to the CS flux) gives an integer 
number of particles per unit cell as required for the ex- 
istence of a gap even though there is no translation 
symmetry breaking in the spin density. 

It is important to check whether the plateau states are 
robust to fluctuations beyond mean-field. Following an 
idea of Fradkin [Q Yang et al. |jl^ showed that Gaussian 
fiuctuations of the CS gauge field are massive provided 
that the TKNN jH] integer describing the quantized Hall 
coefficient of the fermions on the frustrated lattice differs 
from the continuum value of unity. We have computed 
the TKNN integers by following the evolution of gaps 
from the square lattice case where the TKNNs are given 
by a Diophantine equation. None of the plateaus is sup- 
pressed by fluctuations since we find a = —3, —2 and —1 



i and i (see Fig. |). 



respectively for the plateaus at 

At M/ Msat = we have cr = -1 (resp 0) for J / J' < ^(2) 
(resp J /J' > 7(2)). 

A consequence of these non-trivial quantized Hall co- 
efficients for the fermions is that the spin state is chiral 
and exhibits a 'spin quantum Hall effect' |2^. Whether 
this new physics is actually occurring in SrCu2 (603)2 or 



is an unrealistic feature of our mean field approximation 
remains to be seen. One test of our model is the predic- 
tion of a strong plateau in the magnetization beginning 
at GOT, a field which is within reach of modern pulsed 
magnets. 
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